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Abstract — SNS is a MATLAB -based software library written 
to aid in the design and analysis of receiver architectures. It uses 
electrical scattering matrices and noise wave vectors to describe 
receiver architectures of arbitrary topology and complexity. It 
differs from existing freely-available software mainly in that the 
scattering matrices used to describe the receiver and its com- 
ponents are analytic rather than numeric. This allows different 
types of modeling and analysis of receivers to be performed. 

Non-ideal behavior of receiver components can be parameter- 
ized in their scattering matrices. SNS enables the instrument 
designer to then derive analytic expressions for the signal 
and noise at the receiver outputs in terms of parameterized 
component imperfections, and predict their contribution to re- 
ceiver systematic errors precisely. This can drive the receiver 
design process by, for instance, allowing the instrument designer 
to identify which component imperfections contribute most to 
receiver systematic errors, and hence place firm specifications 
on individual components. Using SNS to perform this analysis 
is preferable to traditional Jones matrix-based analysis as it 
includes internal reflections and is able to model noise: two effects 
which Jones matrix analysis is unable to describe. 

SNS can be used to model any receiver in which the com- 
ponents can be described by scattering matrices. Of particular 
interest to the sub-mm and terahertz frequency regime is the 
choice between coherent and direct detection technologies. Steady 
improvements in mm and sub-mm Low Noise Amplifiers (LNAs) 
mean that coherent receivers with LNAs as their first active 
element are becoming increasingly competitive, in terms of 
sensitivity, with bolometer-based receivers at frequencies above 
~ 100 GHz. 

As an example of the utility of SNS, we use it to com- 
pare two polarimeter architectures commonly used to perform 
measurements of the polarized Cosmic Microwave Background: 
differencing polarimeters, an architecture commonly used in 
polarization sensitive bolometer-based polarimeters; and pseudo- 
correlation polarimeters, an architecture commonly used in 
coherent, LNA-based, polarimeters. We parameterize common 
sources of receiver systematic errors in both architectures and 
compare them through their Mueller matrices, which encode 
how well the instruments measure the Stokes parameters of the 
incident radiation. These analytic Mueller matrices are used 
to demonstrate the different sources of systematic errors in 
differencing and correlation polarimeters. 



I. Introduction 

Many fields of astrophysics aim to measure increasingly 
faint signals. For instance, there is great interest at present in 
detecting and characterizing the B-mode [ 1 1 of the polarized 
Cosmic Microwave Background (CMB). The strength of this 
signal is not yet determined by theory, but a strong upper limit 
is 170 nK J2|, a tiny fraction of the CMB total intensity signal 
(~ 2.7 K). 



An instrument built to detect very faint signals will almost 
certainly be heavily affected by systematic errors. It is increas- 
ingly important to be able to model the effects of receiver 
systematic errors on the measured signal, and on the receiver 
sensitivity. We want to be able predict the level of receiver 
systematic errors, show their impact on the gathered data, 
and make quantitative comparisons between different receiver 
architectures. 

Previous analytic and semi-analytic approaches to charac- 
terizing systematic effects in receivers have usually employed 
Jones matrices to describe receiver components and Mueller 
matrices to characterize the effects of receiver systematics on 
the observed Stokes parameters |3|, |4|. The use of Jones ma- 
trices to describe individual receiver components has several 
shortcomings. Only the forward path of the signal through 
the instrument is modeled - internal scattering caused by 
reflections from poorly matched components is not included; 
and Jones matrix modeling is unable to describe component 
noise, and hence receiver sensitivity. Modeling a receiver with 
a full analytic description of the outputs and sensitivity in 
terms of individual component parameters allows us to identify 
which parameters of each component in a receiver are most 
important, and concentrate our efforts on improving them. 

This paper introduces a technique and software for de- 
veloping full analytic descriptions of receiver outputs and 
sensitivities in terms of lab-measureable errors in individual 
components. In this technique components are modeled by 
electrical scattering matrices. When describing a network of 
components with Jones matrices the forward-path cascaded 
response can be obtained through simple matrix manipulation 
and multiplication. The scattering matrix formulation does not 
share this simplicity of calculation: only the case of cascaded 
2-port devices is amenable to a relatively simple analytic 
solution. This paper describes an algorithm for calculating 
the response of arbitrarily connected networks of components. 
We present software which implements this algorithm, and 
apply it to two common polarimeter architectures: differencing 
polarimeters, and pseudo-correlation polarimeters. 

This software allows us to make robust analytic compar- 
isons of receiver architectures. Errors in individual receiver 
components can be parameterized and propagated into the 
description of the receiver performance, e.g. the instrument 
Mueller matrix. We hence have a powerful tool for guiding 
the instrument design process and diagnosing the causes of 
non-ideal instrument behavior. 
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II. Electrical Scattering Matrices 

We model the behavior of individual receiver components, 
and the full receiver, using electrical scattering matrices. 
The electrical scattering matrix (hereafter referred to as the 
scattering matrix) is a representation of a network using the 
ideas of incident, reflected, and transmitted waves. It provides 
a complete description of an Alport network as seen at 
its iV ports. A significant advantage of modeling receiver 
components with scattering matrices is that noise can easily 
be included in the modeling. The noise produced by a device 
is modeled with a noise wave vector; see e.g. 0. 

Consider the arbitrary Alport network shown in Figure la 
We denote the incident wave at port i by V^, the reflected 
wave by , and the noise wave produced by the network at 
that port by Cj. These quantities are related by the scattering 
matrix S and noise wave vector c as follows: 



(1) 



The noise wave voltages Cj of an Alport network are 
complex time-varying random variables characterized by a 
correlation matrix C 
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where the angle brackets indicate time averaging, f indicates 
the conjugate transpose operation, and <g> indicates the outer 
product (or Kronecker product). The diagonal terms of C give 
the noise power deliverable at each port in a 1 Hz bandwidth. 
The off-diagonal terms are correlation products. The noise 
correlation matrix C for a passive network is determined from 
its scattering matrix S by ||6) 



C = fcT(I-SS t ) 



(2) 



where k is Boltzmann's constant, T is the physical temperature 
of the network, and I is the identity matrix. The noise 
correlation matrix for an active network can be determined 
by measurement or modeling. 

III. Solving Arbitrary Networks 
Consider the arbitrarily connected network of AT-port net- 



works shown in Figure lb We need an algorithm to calculate 
the scattering matrix S and noise wave vector c which de- 
scribe the connected network. The algorithm derived here is 
an extension of the algorithm described in j7), with added 
noise wave vector manipulation. Similar algorithms are used 
numerically in SUPERMIX 10. 



First, let us consider the effect of connecting together ports 
k and m of an Alport network described by Equation [T] 
Connecting the ports means that V fc + = and V fe ~ = V£- 
Manipulation of rows k and m of Equation [T] then gives us 
the expressions 
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(4) 



By substituting Equations [3] and |4] into each other we 
can obtain expressions for and VC . Substituting these 
expressions into Equation [T] we can obtain a new expression 
for the reflected wave V~: 
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From Equation |5] we obtain replacement expressions for the 
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Rows and columns k and m are then removed from S and 
rows k and m are removed from c to create the scattering 
matrix and noise vector which describe the new (N — 2)-port 
network. 

A. Algorithm 

To obtain the scattering matrix and noise wave vector which 



describe the arbitrarily connected network shown in Figure lb 
begin by forming the scattering matrix and noise wave vector 
which describe the unconnected network: 
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Fig. 1. (a) An arbitrary TV-port network. The total signal at each port pi 
consists of an incident signal K , a reflected signal V~ , and a noise signal 
d . The reflected signal V~ is a weighted sum of the incident signal at port 
i and transmitted signals from the other ports of the network, the coefficients 
of the sum being the elements of the network scattering matrix S. The noise 
signals are given by the noise wave vector c. (b) An arbitrarily connected 
network of N-port devices. 



We then successively form each connection in the network. 
For each connection, find the rows and columns k and m of 
S and c in Equation [8] which correspond to the ports being 
connected. Use the replacement formulae given by Equations [6] 
and [7] to adjust the S matrix and c vector. Remove rows 
and columns k and m from S, and rows k and m from c. 
Repeat for each remaining connection until we are left with 
the scattering matrix S and noise vector c which describe the 
fully connected network. 

IV. Software Implementation 



We have derived an algorithm in ^III-A for finding the 
scattering matrix and noise wave vector which describe an 
arbitrarily connected network. We need to be able to apply it 
to arbitrary receivers with parameterized scattering matrices 
describing the receiver components and obtain analytic ex- 
pressions for the outputs and noise in terms of the component 
parameters. 

The algorithm must be implemented in a programming 
language with the ability to manipulate symbolic algebraic 




Node 



Fig. 2. Schematic showing the nature of node and TV-port objects and how 
they connect to each other in the software. 



expressions. This programming language must also support 
pointers (or equivalent data structure) to allow the creation 
of a navigateable network. We implemented the algorithm in 
matlabQ which has a powerful and well developed sym- 
bolic algebra toolbox. While it does not have a native pointer 
data type (as of version R2008b), a third party open-source 
pointer librar^ adds this capability. The software package we 
developed to perform this modeling is called SNE 

A. Representing a Network 

A network is represented by nodes and A-port objects, as 
shown in Figure [2] They are both pointer objects. Each A-port 
object contains an array of references to the nodes connected 
to each of its A ports, a scattering matrix S, a noise wave 
vector c, and a variable A, the number of ports of the object. 

Each node object contains two pointers; these refer to the 
A-port objects the node connects to in the "forward" and 
"backward" directions, and which port number the connection 
is made to (Nf p and Nb p respectively). Note that the forward 
and backward directions are completely arbitrary; they are 
merely a helpful concept when trying to visualize the operation 
of the algorithms which act on the network. 

The network is constructed by creating all the node and 
A-port objects using functions called makeNode ( ) and 
makeNport ( ) , assigning scattering matrices and noise wave 
vectors to the A-port objects, and connecting each node to 
its forward and backward A-port objects. A connectNode 
function hides the complexity of assigning references to the 
appropriate array locations and assigning the appropriate port 
numbers to variables. 

Nodes are classified into four types: input, output, central, 
and terminated. We want to calculate the performance of a 
network in terms of the response seen at the outputs due to 
signals presented at the inputs. The central and terminated 
nodes are removed by the network-solving program. 

Once all the objects have been created, assigned matrices 
and vectors, and connected, it suffices to describe the network 
by the four arrays of nodes. Due to the fully connected nature 
of the network representation it is possible to start at any node 
and navigate to any other node by following the appropriate 
links between objects. 

1 http://www.mathworks.com/ 
2 http://code. google, com/p/pointer/ 



3 Download at http://www.astro.caltech.edu/~ogk/SNS/ 
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B. Solving a Network 

The network-solving program accepts four arrays of nodes: 
the inputs, outputs, central nodes, and terminated nodes. It 
returns the scattering matrix and noise wave vector for the 
connected network, where the central and terminated nodes 
have been removed through application of the algorithm 
described in §M-A| 

The first step the software performs is to remove the 
terminated nodes, if there are any. The software assumes 
that the terminations are perfectly matched and at a common 
physical temperature. It modifies the object scattering matrices 
to remove the terminated nodes, and adds the noise terms 
produced by the terminated nodes to the noise wave vectors. 

The network sans terminated nodes is then passed to a 
recursive network shrinking program. This program begins 
with the first central node and applies the algorithm given by 
Equations [6] and [7] to a sub-network consisting of the two N- 
port objects connected to that particular node. A new iV-port 
object is created and assigned the resulting scattering matrix 
and noise wave vector. All the nodes which were connected to 
the now-defunct iV-port objects are reconnected to this new Al- 
port object at the appropriate ports. A new network is formed 
by excluding the central node just considered and the program 
is recursively called on this new network. This continues until 
there are no more central nodes, at which point the scattering 
matrix and noise wave vector of the single remaining Alport 
object are returned. 

Applying the algorithm in this fashion, rather than to the 
entire network at once, means that the size of the matrix S 
in Equation [8] is kept small, speeding up computation. This is 
not an optimum solution to network shrinking, but we have 
found it to be significantly faster than applying the algorithm 
to the full unconnected network for networks of more than a 
few Alport objects. 

The description given above glosses over the significant 
complexity in keeping track of which nodes should be con- 
nected where, and certain configurations of nodes and Alport 
objects which would cause the default implementation of the 
algorithm to fail. The majority of the code is dedicated to 
performing these functions; only a small fraction of the code 
actually carries out the calculations described by the algorithm. 



C. Example Code Listing 

To illustrate the operation of the program consider the 
network shown in Figure [3l It is represented in software by 
connected lists of nodes ana Alport objects, as shown in the 
following code listing: 



% Matrice 
% have be 
% Make the N-po 
PI = makeNport ( 
P4 = makeNport ( 
% Assign scatte 
Pl.S = SI; P2.S 
Pl.c = cl; P2.c 
% Make the node 
Inl = makeNode ( 
0n2 = makeNode ( 
Cn2 = makeNode ( 
Cn4 = makeNode ( 
% Connect nodes 
connectNode (Inl 
connectNode (Cnl 
connectNode (Cn3 
connectNode (Cn5 
connectNode (Onl 
% Build arrays 
inputs = {Inl I 



symbolic algebra 
akeNport ( ) ; 



S2, S3, S4, S5 
viously defined 
rt objects 
; P2 = makeNport (); P3 = : 
; P5 = makeNport ( ) ; 
ing matrices and noise vectors 
= S2; P3.S = S3; P4.S = S4; P5.S = S5; 
= c2; P3.c = c3; P4.c = c4; P5.c = c5; 
which connect the N-port objects 
); In2 = makeNode (); Onl = makeNode (); 
); Inl = makeNode (); Cnl = makeNode (}; 
) ; Cn3 = makeNode O ; 
) ; Cn5 = makeNode () ; 
to N-ports 

, [],1,P1,D; 

PI, 2,P3, 1) ; 
P3, 2,P2, 1) ; 
,P3, 3,P5, 1) ; 
,2, [1,1)1 

; outputs = {Onl 0n2 



ctNode (In2, [ ] , 1,P1, 4) 
ctNode (Cn2,Pl, 3,P4, 1) 
ctNode (Cn4,P4, 2,P2, 4) 
ctNode (Inl, P2, 3, [ ] , 1) 
ctNode (On2,P2, 2, [ ] , 1) 




Fig. 3. An arbitrary network of iV-port devices to illustrate the software. 
Nodes are indicated by open circles, JV-port devices by rectangles. iV-port 
device i is described by scattering matrix S ; and noise wave vector c^. Input 
nodes to the network are In^, central nodes Cn^, terminated nodes Tn^, and 
output nodes On^. 



odes = {Cnl Cn2 Cn3 Cn4 Cn5); tnodes = {Tnl}; 
Pass arrays of nodes to network calculator 

, c] = get ScatteringRecursive (input s , output s, cnodes , tnodes ) ; 

The getScatteringRecursive program performs the 



actions described in { IV-B and returns the scattering matrix S 
and noise wave vector c for the resulting 4-port object. Nodes 
Ini, In2, Oni, and On2 are connected sequentially to ports 1 
to 4 of this object. 

V. POLARIMETRY 



The example presented in the coming section, {VI 



pares two receiver architectures commonly used to measure 
linear polarization at radio to sub-mm wavelengths. This sec- 
tion provides necessary background information by sketching 
a brief summary of polarization. It shows how a receiver may 
be described by a Mueller matrix, and shows how to derive a 
receiver Mueller matrix from the scattering matrix produced 



by the software described in I IV 



A. Brief Summary of Polarization 

An electromagnetic signal is said to be polarized if there is 
some lasting amplitude or phase relation between its orthog- 
onal modes. The coherency vector [9j captures this relation: 

E x (t)E*(ty 
E x {t)E* y {t) 

Ey(t)E*{t) 

E y {t)E* v (t) 
=<Ec3E*) 

Here E is the complex vector of the orthogonal modes E x (t) 
and E y (t) of the signal, (...) indicates time averaging, and <g) 
indicates the outer product. 

If the signal E is acted on by an object described by a Jones 
matrix J, i.e. E ouf = JE, then the new coherency vector is 
given by 



=(J®J*)e 



(9) 



The polarization state of a signal is usually described by 
the Stokes parameters, /, Q, U, and V. I describes the total 
intensity of the signal, Q and U describe the linear polarization 
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state, and V describes the circular polarization state. The 
Stokes vector e s is obtained from the coherency vector e by 



e = 



where T 
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(10) 



(11) 



We see that T is a coordinate transformation of the coherency 
vector to the abstract Stokes frame. 

The Stokes parameters are a convenient and powerful way 
of the describing the state of polarization of an electromagnetic 
signal. From Equation [T0| we have: 

I =(\E x (t)\ 2 ) + (\E y (t)\ 2 ) 
Q={\E x (t)\ 2 )-{\E y (t)\ 2 ) 
U =2^{E x (t)E* y {t)}) 

=(E x (t)E* y (t)) + {E* x (t)E y (t)) 
V =2{%{E x (t)E* y (<)}) 

= - i[{E x (t)E* y (t)) - (E* x (t)E y (t))\ (12) 

a) Mueller Calculus: Suppose that a signal defined by 
the complex electric field vector E and coherency vector e is 
modified by an object described by the Jones matrix J. From 
Equations [9] and 10 we see that the output signal E ont = JE 
will be described by the Stokes vector 

=Me s 

The matrix M = T(J <gi J*)T _1 is called the Mueller matrix. 
It represents the action of the object characterized by Jones 
matrix J in the Stokes vector space. 

Mueller calculus is a matrix method for manipulating Stokes 
vectors. We denote the Mueller matrix elements as 
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Mueller matrices are a convenient means of describing the 
action of an astronomical polarimeter. Of particular interest are 
the Mqi and Mui parameters, which describe the leakage of 
the total intensity / into the measured linear polarization vector 
components. Much of the radio and mm/sub-mm spectrum is 
only slightly linearly polarized, hence non-zero values of Mqi 
and Mui can imply serious contamination of the measured 
linear polarization vector by the total intensity signal. 

B. Deriving Receiver Mueller Matrix 

Say we have calculated the scattering matrix which de- 
scribes the behavior of a receiver. For polarimeters, a natural 
way of expressing the receiver's performance is with a Mueller 
matrix. We need to translate the receiver scattering matrix into 
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Fig. 4. A arbitrary receiver, where orthogonal linear polarizations E x (t) and 
E y (t) are presented at inputs 1 and 2 respectively. D is the output at port 
m. The receiver is described by scattering matrix S. 



a Mueller matrix which describes the action of the receiver on 
the Stokes vector of the incident electromagnetic signal. 

Consider the arbitrary receiver shown in Figure |4] Orthog- 
onal linear polarizations E x (t) and E y (t), representing either 
signals in transmission lines, orthogonal modes in waveguide, 
or orthogonal modes in free-space, are connected to ports 
1 and 2 of the receiver respectively. Receiver output D 
is connected to port m. The receiver is described by the 
scattering matrix S. 

The output E m (t) seen at port m is given by (assuming that 
the connections to ports 3 to N are reflectionless): 

E m (t) =S m iE x (t) + SmzEyit) 

The power contained in the signal E m {t) is then measured. 
At radio wavelengths this is often achieved through the use of 
a square-law detector diode. At mm and sub-mm wavelengths 
a bolometer might be used. The measured power Pr> is given 
by: 

P D =a(E m (t)E m (t)*) 

=a[(|^(i)| 2 )|5 ml | 2 + (|^(t)| 2 )|^ l2 | 2 
+ (E x (t)E* y (t))S ml S; n2 
+ (E* x (t)E y (t))S* ml S m2 \ (13) 

where (...) indicates time averaging, a is a proportionality 
constant dependent on the power detection method, and we 
assume that the instrument scattering matrix parameters are 
constant during the average time period. Now let 

P D =M DI I + M DQ Q + M DU U + M DV V 
=M m {\E x {t)\ 2 + \E y {t)\ 2 ) 
+ M DQ (\E x (t)\ 2 ~ \E y (t)\ 2 ) 
+ M DU (E x (t)E* y (t) + E* x (t)E y (t)) 
- iM DV (E x (t)E* y (t) - E* x (t)E y (t)) (14) 

where we have used the definition of the Stokes parameters 
given in Equation [12| 



By comparing Equations 13 and 14 we can obtain the 
contribution of each Stokes parameter to the power measured 
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at output D: 



Mm 


^{\S m i\ 2 + 
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s m2 \ 2 } 




M DQ 


= f{|5 ml | 2 - 


Sm2\ 2 } 




M DU 








M DV 




~ S ml S m2 } 


(15) 



We can derive the receiver Mueller matrix by applying this 
technique to all the outputs of the receiver. 

VI. POLARIMETER ARCHITECTURE COMPARISON 

Two basic types of architectures are used to measure the po- 
larization of an electromagnetic signal: differencing polarime- 
ters, an architecture commonly used in polarization sensitive 
bolometer-based polarimeters |10|; and pseudo-correlation (or 
correlation) polarimeters, an architecture commonly used in 
coherent, LNA or mixer based, polarimeters such as QUIET 

CD. 

Differencing polarimeters measure the difference in power 
between orthogonal linear modes of the electromagnetic sig- 



nal; see the definition of Q in Equation 12 for inspira- 
tion. Correlation, or pseudo-correlation, architectures measure 
the polarization state by measuring the correlation between 
orthogonal modes. Correlation polarimeters are required to 
preserve the phase of the incident signal. They are hence only 
feasible if coherent (i.e. phase-preserving) amplifiers or mixers 
are available. 

The choice of which architecture to use for a particular 
experiment is often dominated by sensitivity considerations. At 
low frequencies (<~ 60 GHz) the availability of low-noise co- 
herent amplifiers has favored correlation architectures 11121 . At 
higher frequencies, the fundamental quantum limits that ampli- 
fier noise is subject to has favored direct detection technologies 
such as bolometers, and hence differencing polarimeter ar- 
chitectures, for continuum polarimetry experiments. However, 
continuing improvement in coherent amplifier technology at 
high frequencies has pushed their performance closer to the 
quantum limit, e.g. lfl"3l . As coherent amplifier technology 
improves, sensitivity may no longer be the deciding factor 
between technologies, and hence architectures. 

Differencing and correlation architectures measure the po- 
larization information of the incident signal in very different 
ways, and hence suffer from different sources of systematic 
error. A careful analysis of the fundamental strengths and 
weaknesses of each architecture is needed. SNS is well suited 
to perform this analysis. In this section we use SNS to 
derive receiver Mueller matrices for examples of these two 
polarimeter architectures. This analysis highlights the different 
sources of systematic error in these architectures. 

A. Differencing Polarimeters 

An example of a differencing polarimeter architecture is 
shown in Figure [5] (right). The powers in orthogonal linear 
modes Di and D 2 are detected and differenced to obtain one 
of the linear polarization parameters. Differencing polarime- 
ters have a much simpler architecture than pseudo-correlation 
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Fig. 5. Examples of two commonly used polarimeter architectures, (left) A 
pseudo-correlation polarimeter. (right) A differencing polarimeter. Orthogonal 
waveguide modes E x and E y are circularized or rotated, and extracted from 
waveguide by an orthomode transducer (OMT). In the pseudo-correlation 
architecture the signals are further processed. The powers in signals Di to 
D4 are detected and processed to obtain the Stokes parameters, as explained 
in the text. 



polarimeters. However, they measure only a single linear 
Stokes parameter; a duplicate receiver oriented at 45° to the 
first is needed to measure the second linear Stokes parameter. 



B. Pseudo-correlation Polarimeters 

An example of a pseudo-correlation polarimeter architecture 
is shown in Figure [5] (left). The Stokes parameters of linear 
polarization in a circular basis are given by: 



Ei(t) 



[E X (t) - lEy(i)\ 



1 



E r (t)=—[E x (t)+iE y (t)] 

Q=2(n{E l (t)E* r (t)}) 
U = -2(Z{E l (t)E* r (t)}) 



(16) 



Pseudo-correlation polarimeters measure the linear Stokes 
vector by correlating circular polarization signals Ei and E r 
with (Z?3 and D4) and without (Di and D 2 ) a 90° phase shift. 
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The powers in signals D\ to D4 are detected and combined 
to obtain the linear Stokes parameters. 

While the architecture of a pseudo-correlation polarimeter is 
more complex than that of a differencing polarimeter, it does 
provide some significant advantages. For instance, both linear 
polarization parameters can be measured with a single optical 
assembly, providing twice the information for the same focal 
plane area occupied. 

C. Parameterized Scattering Matrices 

Some of the components in the receivers shown in Figure [5] 
have parameterized scattering matrices. While it is possible to 
describe every component in a receiver with a suitable parame- 
terized scattering matrix, the resulting analytic expressions for 
the outputs soon become too complicated to be useful when 
written down. 

In this analysis the components are assumed to be perfectly 
matched, i.e. the diagonal elements of the parameterized 
scattering matrices are zero. 

1) Circularizes The circular phase shifter translates or- 
thogonal linear polarizations into orthogonal circular polar- 
izations. It introduces a 90° phase shift into one orthogonal 
linear mode, and is oriented at 45° to the OMT linear axis. 

A possible parameterization of the circularizer's scattering 
matrix as shown in Figure [5] is: 

ro 1 1 
1 
1 

_ e i(f+9o) e »(§+0e) 



Lc 
V2 





. e »(f+fc) 







Here L? c is the insertion loss of the circularizer, and 8 C is the 
error in the 90° phase shift. The circularizer is assumed to be 
otherwise perfect. 

2) Faraday Rotator: The Faraday rotator shown in the 
differencing polarimeter in Figure |5]is a component sometimes 
used in PSBs [14|. It modulates the measured polarization 
signal by introducing a variable rotation to the plane of linear 
polarization of the incident signal. A similar effect can be 
achieved with a rotating birefringent half-waveplate, or a wire 
grid. 

The scattering matrix for the rotator considered here is: 





cos(9 ps ) 
sm(9 ps ) 
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— sinfi 



,) 
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ps 
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Here 28 ps is the time-dependent linear plane rotation intro- 
duced by the Faraday rotator. 

3) OMT: The OMT extracts orthogonal linear modes from 
the waveguide. The scattering matrix parameterization consid- 
ered here is: 



'OA/T 



Here D x and D y measure the insertion loss for each orthogo- 
nal mode, while d xy and d yx measure the leakage of one mode 
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dyx 





D x 








dxy 
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dxy 
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into the other. From conservation of energy considerations in 
a passive component we have the constraints \D. 



\D\, \d 



xy\ 



*yx\ 



% and \D\ 



A = \Dy\ = 

Llmt> where 



L 2 omt is the insertion loss of the OMT. The parameters may 
have arbitrary phase. 

4) LNAs and Phase Switch: The scattering matrices for the 
LNAs in the pseudo-correlation polarimeter are given by: 



&l,r — 





Gl,r 



Here Gl and Gr are the complex voltage gains of the left 
and right circular polarization amplifiers respectively. 

The phase switch in the pseudo-correlation polarimeter 
modulates the phase of one of the signal arms relative to the 
other. Its scattering matrix is given by: 





A8„ 







Here 9 ps is the time-dependent phase shift introduced by the 
phase switch, usually shifting between 0° and 180°. 

D. Polarimeter Mueller Matrix Elements 

We now build a connected model of each polarimeter in 



SNS using the matrix parameterizations given in £ VI-C We 



use Equation 15 to obtain expressions for the powers measured 
at each receiver output in terms of the incident signal Stokes 
parameters. 

For the ideal pseudo-correlation polarimeter, i.e. where all 
the components are perfect, the outputs are: 
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P D1 
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Pd 3 
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P Di 







(17) 



Here we have assumed that 6 ps = and a = 1. To measure 
the Stokes parameters we take Q m = Pp A — Pd 3 , U m = 
Pd 2 - Pdi, and I m = \{P Dl + P D2 + P Ds + P Di ). 
For the ideal differencing polarimeter the outputs are: 



P Dl =^I+^cos(29 ps )Q 
Pd 2 =\l - l -cos{26 ps )Q 



sm(26 ps )U 
sin(20 ps )C7 



(18) 



Which linear Stokes parameter we measure depends on the 
plane rotation introduced by the Faraday rotator, and is given 
by L rn = P Dl - P D . 2 = cos(26 ps )Q - sm{20 ps )U. The 
measured total intensity is I m = Pj) 1 + Pd 2 - 

The Mueller matrix parameters of particular interest in CMB 
polarization studies are: M/j, Mqq, and Muu> me diagonal 
elements of the Mueller matrix; Mqi and Mui, the leakage 
of the total intensity signal into the (generally) small linear 
polarization signal; and Mqjj and Mjjq, which measure the 
rotation of the linear polarization vector by the receiver. 
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These Mueller matrix parameters for the pseudo-correlation 
polarimeter are: 



Mr, = - 



[I 



G L 



\Gi 



1 



M QQ =L 2 c \G l G r \ [|D! 2 cos(6» 3 ) + jdj 2 cos(0 4 )] 

£ 2 |GlG h |[cos(0 c ){|Z)| 2 cos(6» 3 ) — ja!| 2 cos(6 l 4)} 
- sin(6» c )|Dd|{sm(6'i) + sin(0 2 )}l 
L 2 c \G l Gr\\ Dd\ [ cos(6» 2 ) + cos(0i )] 
L 2 |G L G fl ||Dd|[sin(e 1 )-sin(e 2 )] 
M QV =L 2 |G£GH|[cos(e c ){|Di 2 sin(6)3)- |d| 2 sin(r? 4 )} 
+ sin(0 c )|_Dd|{cos(6> 2 ) - cos(0i)}] 

L 2 |G£Gfl| [|D| 2 sin(0 3 ) + |d| 2 sin(0 4 )] (19) 



Muu =1 
Mqi 



Muq 
where 9\ 

e 3 



9d v - 0d xy - (8g l - 0g r + 6 P s) 
- dd yx + (0 Gl - 8 Gr + d ps ) 

9 Dx -e Dy + {e Gh -o GR +o ps ) 

9d xy - d + (0 Gl - 6 Gr + 6 ps ) 



Here X = \X\e x . We have implicitly assumed that the 
responses of all the power detectors are equal and stable. 

To keep the comparison between the architectures reason- 
able, we need to include varying power detection sensitivity in 
the differencing polarimeter. Let the power detection propor- 
tionality constants (see Equation [D) be u\ and for outputs 
1 and 2 respectively. We also neea to decide on the rotation 
angle of the Faraday rotator to specify which linear Stokes 
parameter we actually measure. Let 8 ps = ±45° (i.e. Stokes 
U). We then obtain the Mueller matrix parameters: 



M vu 

Mui 
Muq 



Lomt \ | 1 
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sin(26> ps ) r ] r 
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— Ct2 cos(6d 



\ol2 — Qij ^ 



Ql COS! tfn — U dx 



D \ - \ d 
— Ct2 cos(6d 



cos(29 ps )\Dd\ 
- sm(26 ps )\Dd\ 



(20) 



1 ) Discussion: One of the greatest sources of systematic 
error in polarimeters is leakage of the total intensity signal into 
the measured linear polarization amplitude, P = \/Q 2 + U 2 . 
The fractional contribution to P from total intensity leakage, 
APj = Mpj/Mu, is given by: 



M 2 
m QI 



Assume that we have two differencing polarimeters oriented 
such that they measure Q and U respectively, identical ex- 
cept for their values of a. The "Q" polarimeter has values 
ai, ct2, while the "U" polarimeter has a^, a^. For the pseudo- 
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Fig. 6. Comparison of the effect of imbalance on the fractional total intensity 
to polarization leakage for pseudo-correlation and differencing polarimeter 
architectures. For the pseudo-correlation architecture, A = A|G| 2 . For the 
differencing architecture A = Aa. \d\ 2 is the OMT cross polarization. 



correlation and differencing polarimeters we then have: 

2V2\G L G R \ 



APf = 



APf 



G L \ 2 + \G R , 
J(ai- a 2 ) 2 + (a 3 



d\V(l-\d\ 2 )(l 
ai) 2 



cos 



0) 



a.2 + 013 + oti 

As a simplification, assume that Aa = (ai — 02) /a2 = («3 — 
a 4 )/a 4 , and let A|G| 2 = ^feif ^ ■ Take the worst-case 
phase scenario, where cos(#i +82) = 1. We now have: 

y/l + A|G| 2 
2 + AIGI 2 



APf C =4|d|VTHdp 



APf 



V2Ac 



(21) 

2 + Aa (22) 
Equations 21 and 22 are plotted in Figure [6] Several 
attributes are noteworthy: APj is independent of the OMT 
cross polarization \d\ 2 for the differencing polarimeter, but 
is heavily dependent on the power sensitivity imbalance Aa; 
APj is almost independent of gain imbalance A|G| 2 for the 
pseudo-correlation polarimeter, but is dependent on the OMT 
cross polarization. 

Figure [6] clearly illustrates the difference between the po- 
larimeter architectures in terms of total intensity to polarization 
leakage. Correlation polarimeters are very insensitive to what 
is generally the most unstable parameter in a coherent receiver: 
fluctuating LNA gain. They are moderately sensitive to OMT 
cross polarization |g?| 2 . The comparatively high sensitivity of 
differencing polarimeters to power detection imbalance can be 
reduced if Aa is stable and well-known; the data can then be 
corrected and the leakage of / into P reduced. 

E. Pseudo-Correlation Polarimeter Noise Temperature 

A very powerful benefit of using scattering matrices to 
model receivers is the ability to perform noise analysis. We 
specify parameterized noise wave vectors for the pre-LNA 
components and ignore any noise produced by the components 
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"down stream" of the LNAs, as their contribution will be 
negligible if the LNA gain is high. 

If the noise wave vector of the pseudo-correlation polarime- 
ter is given by c, then the noise power measured at the output 
Di in a 1 Hz bandwidth is given by Pj = a(cjC*), where Cj 
is the noise wave vector element corresponding to output Di. 

Suppose that component k of M total components has N 
ports, and is specified by the scattering matrix S fc and the 
noise wave vector c k . Cj is given by: 



M 

fe=i 



c k , where c k = 



N 

£■ 

i=l 



Noise waves from different devices are not correlated: 

-iidj 1 )*) = for k 7^ m. So, Pj is given by: 



A I 



k=l 

where P 4 fe =^^(c fc .(b fe ®(b fc )t; 



Here C k is the noise correlation matrix for component k, ■ 
is the matrix dot product, b fc is the vector [b k . . . &^-] T , <E> is 
the outer product, and t is the hermitian conjugate. The Yl 
indicates a sum over all the matrix elements. 

The noise correlation matrices for the passive pre-LNA 
components are obtained using Equation [2] To get the noise 
correlation matrices for the LNAs we make two simplifications 
to the HEMT noise correlation matrix model in [5|: first, 
the off-diagonal terms of an LNA noise correlation matrix 
((cidjj) and (cjcg)) are much smaller than the diagonal terms 
so we take them to be zero. Second: { | C2 1 2 } — |<S2i| 2 (|ci| 2 ) — 
k\S2i\ 2 Ts, where Tn is the amplifier noise temperature and k 
is Boltzmann's constant. 



Receiver noise temperature is referenced to the input. We 
consider the receiver temperature at output Di to be 
the temperature of a thermal source seen equally at each 
input which produces the same total output noise power in 
a noiseless receiver: 



--a(jT\S k 



.Ti = 



Ti 



P- 



(23) 



Here Aj n is the number of inputs to the receiver, and S is the 
receiver scattering matrix. 

Applying this technique to the pseudo-correlation polarime- 
ter we derive the receiver temperatures for the outputs D\ to 




Fig. 7. Plots of Ti/Tjv against various parameters. We see that the receiver 
is most sensitive to the cross polarization, |d| 2 , and is negligibly sensitive to 
the other parameters. We assume Tjy = Tp = 15 K. See text for details. 
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(24) 



where T c =T V 



1 



G=\G L \ 2 + \G 2 R 



E=(\G L \ 2 + \G 2 R \)Ll mt 

F =2\G L G R \\Dd\ [cos(0i) + cos(6> 2 ) 

H =2\G L G R \\Dd\ [sin(^) - sin(c? 2 ) 



Here we have assumed that the amplifiers have noise temper- 
ature Tn, and that the circularize! - and OMT are at a physical 



temperature of T p . 6\ and 6*2 are given in Equation 19 

1) Discussion: Which parameters in the multiparameter 
expression for T% in Equation [24] have the greatest impact on 
the receiver temperature? The most obvious are L 2 C and L 2 mt , 
the insertion losses of the circularizer and OMT respectively. 
Setting those aside, how do the other parameters affect the 
receiver temperature? 

The minimum value that T\ can have is Tjy. Our sensitivity 
impact metric is then T X /T N > 1. We fix L 2 = —0.1 dB and 
L 2 omt = -0.2 dB, and set \G R \ = 1, \G L \ = y/l + A|G| 2 , 
9q l = 0, and 9d x =0 (only relative phases matter here). This 
leaves us with six parameters: \d\ 2 , A|G| 2 , 0q r , &d , Qd^ , 
and 9d y:c - 

We generate random sets of physically realistic values for 
these parameters, and evaluate the metric Ti/Tjv for each set. 
We plot Ti/T/v against each parameter under consideration 
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in Figure [7] It is immediately clear that the most important 
parameter in this set is \d\ 2 . Discarding the least important 
parameters we now have: 



'LI Ll mt -2\d\^/L- 
T P (1-L 2 C L, 



' 2 T 2 ] 
omt J 



2 

omt 



L 2 c L 2 omt 



since \d\ 2 < L 2 Q 



(25) 



This simplified expression for T\ is exactly what we would 
derive using a conventional noise temperature analysis, indi- 
cating that the software has calculated the noise temperature 
correctly. 

VII. Conclusions 

We have presented SNS, a MATLAB-based software library 
written to aid in the design and analysis of receiver archi- 
tectures. It uses electrical scattering matrices and noise wave 
vectors to describe receiver architectures of arbitrary topology 
and complexity. 

We use SNS to compare two polarimeter architectures com- 
monly used to perform measurements of the polarized CMB: 
differencing polarimeters, an architecture commonly used in 
PSB-based polarimeters; and pseudo-correlation polarimeters, 
an architecture commonly used in coherent polarimeters. This 
analysis highlights the differing sources of systematic error in 
these architectures: I to P leakage in pseudo-correlation po- 
larimeters is almost immune to gain imbalance, but sensitive to 
OMT cross polarization; while I to P leakage in differencing 
polarimeters is immune to OMT cross polarization, but very 
sensitive to power detection imbalance. 

We show how SNS can be used to calculate analytical 
expressions for the receiver noise temperature of arbitrary 
receivers. Analytic expressions for the receiver temperature of 
a pseudo-correlation polarimeter are derived, and are found to 
be consistent with those obtained from conventional receiver 
temperature calculations. 
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